Short-time critical dynamics of the three-dimensional systems 
with long-range correlated disorder 



Vladimir V. Prudnikov 1 '*, Pavel V. Prudnikov 1 , Bo Zheng 2 , 
Sergei V. Dorofeev 1 and Vyacheslav Yu. Kolesnikov 1 , 
1 Dept. of Theoretical Physics, Omsk State University, Omsk 644077, Russia 
^Physics Dept., Zhejiang University, Hangzhou 310027, PR of China 

Monte Carlo simulations of the short-time dynamic behavior are reported for three- 
dimensional Ising and XY models with long-range correlated disorder at criticality, in 
the case corresponding to linear defects. The static and dynamic critical exponents 
are determined for systems starting separately from ordered and disordered initial 
states. The obtained values of the exponents are in a good agreement with results of 
the field-theoretic description of the critical behavior of these models in the two-loop 
approximation and with our results of Monte Carlo simulations of three-dimensional 
Ising model in equilibrium state. 

I. INTRODUCTION 

The investigation of critical behavior of disordered systems remains one of the main 
problems in condensed matter physics and excites a great interest, because all real solids 
contain structural defects. The structural disorder breaks the translational symmetry of 
the crystal and thus greatly complicates the theoretical description of the material. The 
influence of disorder is particularly important near critical point where behavior of a system 
is characterized by anomalous large response on any even weak perturbation. In most 
investigations consideration has been restricted to the case of point-like uncorrelated defects 
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However, the non-idealities of structure cannot be modeled by simple uncorrelated 
defects only. Solids often contain defects of a more complex structure: linear dislocations, 
planar grain boundaries, clusters of point-like defects, and so on. 

Different models of structural disorder have arisen as an attempt to describe such com- 
plicated defects. In this paper we concentrate on model of Weinrib and Halperin (WH) 
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with the so-called long-range correlated disorder when pair correlation function for point-like 
defects g(x — y) falls off with distance as a power law g(x — y) ~ |x — y| _a . Weinrib and 
Halperin showed that for a > d long-range correlations are irrelevant and the usual short- 
range Harris criterion 2 — dv§ = «o > of the effect of point-like uncorrelated defects is 
realized, where d is the spatial dimension, and vq and cto are the correlation-length and the 
specific-heat exponents of the pure system. For a < d the extended criterion 2 — au > 
of the effect of disorder on the critical behavior was established. As a result, a wider class 
of disordered systems, not only the three-dimensional diluted Ising model with point-like 
uncorrelated defects, can be characterized by a new type of critical behavior. So, for a < d 
a new long-range (LR) disorder stable fixed point (FP) of the renormalization group recur- 
sion relations for systems with a number of components of the order parameter m > 2 was 
discovered. The critical exponents were calculated in the one-loop approximation using a 
double expansion in e = 4 — d 1 and 5 = 4 — a <C 1 . The correlation-length exponent was 
evaluated in this linear approximation as v = 2/ a and it was argued that this scaling relation 
is exact and also holds in higher order approximation. In the case m = 1 the accidental 
degeneracy of the recursion relations in the one-loop approximation did not permit to find 
LR disorder stable FP. Korzhenevskii et al [4] proved the existence of the LR disorder stable 
FP for the one-component WH model and also found characteristics of this type of critical 
behavior. 

Ballesteros and Parisi [5] have studied by Monte Carlo means the critical behavior in 
equilibrium of the 3D site-diluted Ising model with LR spatially correlated disorder, in the 
a = 2 case corresponding to linear defects. They have computed the critical exponents of 
these systems with the use of the finite-size scaling techniques and found that a v value is 
compatible with the analytical predictions v — 2/ a. 

However, numerous investigations of pure and disordered systems performed with the use 
of the field-theoretic approach show that the predictions made in the one-loop approxima- 
tion, especially based on the e - expansion, can differ strongly from the real critical behavior 
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Therefore, the results for WH model with LR correlated de: 
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13], where a renor- 



121 ] was questioned in our paper 
malization analysis of scaling functions was carried out directly for the 3D systems in the 
two- loop approximation with the values of a in the range 2 < a < 3, and the FPs corre- 
sponding to stability of various types of critical behavior were identified. The static and 



dynamic critical exponents in the two-loop approximation were calculated with the use of 



the Pade-Borel summation technique. The results obtained in [13j essentially differ from the 
results evaluated by a double e, 5 - expansion. The comparison of calculated the exponent v 
values and ratio 2/ a showed the violation of the relation v — 2/ a, supposed in j^j] as exact. 

The models with LR-correlated quenched defects have both theoretical interest due to 
the possibility of predicting new types of critical behavior in disordered systems and exper- 
imental interest due to the possibility of realizing LR-correlated defects in the orientational 
glasses [ijj], polymers [if] ], and disordered solids containing fractal-like defects Q] or dislo- 
cations near the sample surface [if]]. 

To shed light on the reason of discrepancy between the results of Monte Carlo simulation 



of the 3D Ising model with LR-correlated disorder 



our renormalization group description of this model 
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in the a = 2 case, and the results of 
we have computed by the short- 



time dynamics method [l3, Q the static and dynamic critical exponents for site-diluted 3D 
Ising and XY models with the linear defects of random orientation in a sample. 

In the following section, we introduce a site-diluted 3D Ising model with the linear defects 
and scaling relations for the short-time critical dynamics. In Section IHIl we give results of 
critical temperature determination for 3D Ising model with the linear defects for case with 
spin concentration p = 0.8. We analyze the critical short-time dynamics in Ising systems 
starting separately from ordered and disordered initial states. Critical exponents obtained 
under these two conditions with the use of the corrections to scaling are compared. Also, 
in Section IHIl the results of measurements of the critical characteristics for 3D Ising model 
in equilibrium state are presented in comparison with results of the short-time dynamics 
method. In Section IIVI the results of Monte Carlo studies of critical behavior of 3D XY- 
model with linear defects for the same spin concentration p = 0.8 are considered. Our main 
conclusions are discussed in Section |Vj 

II. MODEL AND METHODS 

We have considered the following Ising model Hamiltonian defined in a cubic lattice of 
linear size L with periodic boundary conditions: 

H = -J ^PiPjSiSj, (1) 
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where the sum is extended to the nearest neighbors, Si = ±1 are the usual spin variables, 
and the Pi are quenched random variables (pi = 1, when the site i is occupied by spin, and 
Pi = 0, when the site is empty), with LR spatial correlation. An actual pi set will be called 
a sample from now on. We have studied the next way to introduce the correlation between 
the pi variables for WH model with a = 2, corresponding to linear defects. We start with a 
filled cubic lattice and remove lines of spins until we get the fixed spin concentration p in the 
sample. We remove lines along the coordinate axes only to preserve the lattice symmetries 
and equalize the probability of removal for all the lattice points. This model was referred 
in [s| as the model with non-Gaussian distribution noise and characterized by the isotropic 
impurity-impurity pair correlation function decays for large r as g(r) ~ 1/r 2 . In contrast to 
j^l we put a condition of linear defects disjointness on their distribution in a sample, whereas 
in (ll the possibility of linear defects intersection is not discarded. The physical grounds for 
this condition are connected with fact that in real materials dislocations as linear defects 
are distributed uniformly in macroscopic sample with probability of their intersection close 
to zero. The condition of linear defects disjointness corresponds to WH model since the 
intersection of linear defects being taken into consideration results in additional vertices of 
interaction which are absent in the effective Hamiltonian of WH model. 

In this paper we have investigated systems with the spin concentration p = 0.8. We have 
considered the cubic lattices with linear sizes L from 16 to 128. The Metropolis algorithm 
has been used in simulations. We consider only the dynamic evolution of systems described 
by the model A in the classification of Hohenberg and Halperin • The Metropolis Monte 
Carlo scheme of simulation with the dynamics of a single-spin flips reflects the dynamics 
of model A and enables us to compare the obtained dynamical critical exponent z with 
the results of our renormalization group description of critical dynamics of this model [3] 
having LR-disorder. 

A lot of results have been recently obtained concerning the critical dynamical behavior of 
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statistical models |17l . 1181 ] in the macroscopic short-time regime. This kind of investigation 
was motivated by analytical and numerical results contained in the papers of Janssen et 
al (2^] and Huse 2l|. Important is that extra critical exponents should be introduced to 
describe the dependence of the scaling behavior for thermodynamic and correlation functions 
on the initial conditions. According to the argument of Janssen, Schaub and Schmittman 
(20I ] obtained with the renormalization group method, one may expect a generalized scaling 



relation for the k-th moment of the magnetization 

M (fc) (t, r, L, mo) = b~ k P /u M (k) (b- z t, b 1/u r, b^L, b x °m ) , (2) 

is realized after a time scale t mic which is large enough in microscopic sense but still very 
small in macroscopic sense. In /3, v are the well-known static critical exponents and z is 
the dynamic exponent, while the new independent exponent xo is the scaling dimension of 
the initial magnetization mo, r = (T — T c )/T c is the reduced temperature. 

Since the system is in the early stage of the evolution the correlation length is still small 
and finite size problems are nearly absent. Therefore we generally consider L large enough 
and skip this argument. We now choose the scaling factor b = t 1 ^ so that the main t- 
dependence on the right is cancelled. Expanding the scaling form §2§ for k = 1 with respect 
to the small quantity t Xo ' z mo, one obtains 

M(t, r, m ) ~ m t 6 F(t l/uz r, t Xo/z m ) = m t e (l + at 1,vz r) + 0(r 2 , m 2 ), (3) 

where 9 = {x$—f3 jv)jz has been introduced. For r = and small enough t and mo the scaling 
dependence for magnetization (jBJ) takes the form M(t) ~ mot 9 . For almost all statistical 



systems studied up to now [13, ll8|, |22j], the exponent 9 is positive, i.e., the magnetization 
undergoes surprisingly a critical initial increase. The time scale of this increase is to ~ 
m - z / x o_ However, in the limit of m the time scale goes to infinity. Hence the initial 
condition can leave its trace even in the long-time regime. 

If r 7^ 0, the power law behavior is modified by the scaling function F{t l l uz T) with 
corrections to the simple power law, which will be depended on the sign of r. Therefore, 
simulation of the system for temperatures near the critical point allows to obtain the time 
dependent magnetization with non-perfect power behavior, and the critical temperature T c 
can be determined by interpolation. 

Other two interesting observables in short-time dynamics are the second moment of mag- 
netization M^it) and the auto-correlation function 

A{t) = ^ d (^S l {t)SM^J- (4) 

As the spatial correlation length in the beginning of the time evolution is small, for a finite 
system of dimension d with lattice size L the second moment M^ 2 \t, L) ~ L~ d . Combining 



this with the result of the scaling form (121) for r = and b = t 1 ^, one obtains 

M (2) (t) _ t -2 P /uz M (2) ( 1)f -l/* L ) „ t C 2> C2 = f d _ I (5) 

lermore, careful scaling analysis shows that the auto-correlation also decays by a power 
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A(t)~t~ c «, c a = - z -9. (6) 

Thus, the investigation of the short-time evolution of system from a high-temperature initial 
state with m = allows to determine the dynamic exponent z, the ratio of static exponents 
/3/u and a new independent critical exponent 9. 

Till now a completely disordered initial state has been considered as starting point, i.e., 
a state of very high temperature. The question arises how a completely ordered initial state 
evolves, when heated up suddenly to the critical temperature. In the scaling form (j2J) one 
can skip besides L, also the argument m = 1 

M (fc) (t, r) = r w "M w (b- z t, b 1/u r) . (7) 

The system is simulated numerically by starting with a completely ordered state, whose 
evaluation is measured at or near the critical temperature. The quantities measured are 
M(t), M®(t). With b = t l l z one avoids the main t-dependence in M^ k \t), and for k = 1 
one has 

M(t, r) = rP /vz M(l, t 1/uz r) = t~ p/vz (1 + at 1/uz T + 0(r 2 )) . (8) 

For r = the magnetization decays by a power law M(t) ~ t- p f vz . If r ^ 0, the power 
law behavior is modified by the scaling function M{l ) t 1 ^ uz T). From this fact, the critical 
temperature T c and the critical exponent /3/uz can be determined. 

We must note, that the short-time dynamic method in part of critical evolution de- 
scription of system starting from the ordered initial state is essentially the same as the 



non-equilibrium relaxation method proposed by N. Ito in [24( for critical behavior study of 
three-dimensional pure Ising model. At present, this method was extended by N. Ito to 
non-equilibrium relaxation study of Ising spin glass models 25[], Kosterlitz-Thouless phase 
transition 26|] and fully frustrated XY models in two dimension 27]. 

From scaling form (jSJ) the power law of time dependence for the logarifmic derivative of 
the magnetization can be obtained in the next form 

9>M(f,r)| T=0 ~^, (9) 



which allows to determine the ratio \jvz. On basis of the magnetization and its second 
moment the time dependent Binder cumulant 

A/f( 2 ) 

Ut) = w - 1 - t d/z (10) 

is defined. From its slope one can directly measure the dynamic exponent z. Consequently, 
from an investigation of the system relaxation from ordered initial state with m = 1 the 
dynamic exponent z and the static exponents (5 and v can be determined and their values 
can be compared with results of simulation of system behavior from disordered initial state 
with m Q = 0. 



III. MEASUREMENTS OF THE CRITICAL TEMPERATURE AND CRITICAL 

EXPONENTS FOR 3D ISING MODEL 

We have performed simulations on three-dimensional cubic lattices with linear sizes L 
from 16 to 128, starting either from an ordered state or from a high-temperature state with 
zero or small initial magnetization. We would like to mention that measurements starting 
from from a completely ordered state with the spins oriented in the same direction (m = 1) 
are more favorable, since they are much less affected by fluctuations, because the quantities 
measured are rather big in contrast to those from a random start with m = 1. Therefore, for 
careful determination of the critical temperature and critical exponents for 3D Ising model 
with linear defects we begin to investigate the relaxation of this model from a completely 
ordered initial state. 



A. Evolution from an ordered state 



Initial configurations for systems with the spin concentration p = 0.8 and randomly 
distributed quenched linear defects were generated numerically. Starting from those initial 
configurations, the system was updated with Metropolis algorithm. Simulation have been 
performed up to t = 1000. We measured the time evolution of the magnetization 
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FIG. 1: Time evolution of the magnetization M(t) for L = 128 and for different values of the 
temperature T. 



and the second moment 
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which also allow to calculate the time dependent Binder cumulant U^it) ffTUj) . The angle 
brackets in fill I) and (I12p denote the statistical averaging and the square brackets are for 
averaging over the different impurity configurations. 

In Fig. [1] the magnetization M(t) for samples with linear size L = 128 at T = 3.919, 
3.925, 3.930, 3.935 and 3.940 is plotted in log-log scale. The resulting curves in Fig. CD have 
been obtained by averaging over 3000 samples with different linear defects configurations. 
We have determined the critical temperature T c = 3.930(2) from best fitting of these curves 
by power law. 

The critical temperature determined in 5] for the same system with spin concentration 
p = 0.8 in the non-Gaussian case is T c = 3.8891(2). This difference of the critical tempera- 
ture values shows that different principles of distribution of linear defects are the reason of 
discrepancy between the results obtained in by Monte Carlo simulation of the 3D Ising 
model with LR-correlated disorder, and results in renormalization group description of this 



model [13j. 

In order to check-up the critical temperature value independently, we have carried out in 
equilibrium the calculation of Binder cumulant U4, defined as 
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and the correlation length 28] 
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where (xi t j,X2j,xsj) are coordinates of j-th site of lattice. 
The cumulant U±(L,T) has a scaling form 

U 4 (L,T) = u(L 1 ^(T-T c )). 



(18) 



The scaling dependence of the cumulant makes it possible to determine the critical tem- 
perature T c from the coordinate of the points of intersections of the curves specifying the 
temperature dependence U^L, T) for different L. In Fig. [2K the computed curves of U^L, T) 
are presented for lattices with sizes L from 16 to 128. As a result it was determined that the 
critical temperature is T c = 3.9275(5). In this case for simulations we have used the Wolff 
single-cluster algorithm with elementary MCS step as 5 cluster flips. We discard 10000 MCS 
for equilibration and then measure after every MCS with averaging over 100000 MCS. The 
results have been averaged over 15000 different samples for lattices with sizes L = 16, 32 
and over 10000 samples for lattices with sizes L = 64, 128. 

The crossing of £/L was introduced as a convenient method for calculating of T c in 



29l | . In Fig. [2b the computed curves of temperature dependence of ratio £/L are presented 
for lattices with the same sizes, the coordinate of the points of intersections of which also 
gives the critical temperature T c = 3.9281(1). This value of T c we selected as the best for 
subsequent investigations of the Ising model. 





FIG. 3: Time evolution of logarithmic derivative of the magnetization d T lnM(t, t)| ( a ) an d 
the cumulant Uz{t) (b) for L = 128 at the critical temperature T c = 3.9281. 

Also, we have determined the temperature of intersection of the curves specifying the 
temperature dependence cumulants U^L, T) for L = 16 and L = 32 with the use of linear de- 
fects distribution in samples as in with the possibility of their intersection. Computation 
gives T C (L) = 3.8884(6) in this case which corresponds to the results in [5(] but differs from 
T C (L) = 3.9185(5) obtained with the use of condition of linear defects disjointness for lattices 
with the same sizes. Turning back to short-time dynamics method, we note that the expo- 
nent l/i/z can be determined from relation © if we differentiate lnM(t, r) with respect to r. 



The dynamic exponent z can be determined from 
analysis of time dependent Binder cumulant U^if) 
(110p for t = 0. In Fig.[3]the logarithmic derivative 
of the magnetization d T lnM(t, t)| t=0 with respect 
to t (Fig. [3h) and the cumulant ^(t) (Fig. [3b) for 
samples with linear size L = 128 at T c = 3.9281 
are plotted in log-log scale. The d T lnM(t,r)\ T=0 
have been obtained from a quadratic interpolation 
between the three curves of time evolution of the 
magnetization for the temperatures T = 3.9250, 
3.9281, 3.9310 and taken at the critical tempera- 
ture T c = 3.9281. The resulting curves have been 
obtained by averaging over 3000 samples. 

We have analyzed the time dependence of the 
cumulant U^if) and clarified that in the time in- 
terval [50,150] the Z72 (t) is best fitted by power 
law with the dynamic exponent z ~ 2.02, corre- 
sponding to the pure Ising model [30|, [3l|, and 



the linear defects are developed for t > 400 MCS 
only. An analysis of the U%(t) slope measured 
in the interval [500,900] shows that the exponent 
d/z = 1.173(12) which gives z = 2.558(26). We 
have taken into account these dynamic crossover 
effects for analysis of the time dependence of mag- 
netization and its derivative. So, the slope of 
magnetization and its derivative over the interval 
[450,900] provides the exponets fi/vz = 0.230(2) 
and 1/vz = 0.517(12) which give u = 0.746(19) 
and = 0.445(10). 

In the next stage, we have considered the cor- 
rections to the scaling in order to obtain accurate 
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FIG. 4: Dependence of the mean square 
errors a of the fits for the magnetization 
(a), logarifmic derivative of the magneti- 
zation (b) , and the cumulant (c) as a func- 
tion of the exponents 0/vz, 1/vz, and d/z 
for uj = 0.8. 



values of the critical exponents. We have applied the following expression for the observ- 



TABLE I: Values of the exponents 0/vz, 1/uz, d/z, and minimal values of the mean square errors 
a in fits for different values of the exponent uj 



uj (3/vz a 1/vz a d/z a 
0.7 0.2112 0.0100 0.556 0.0053 1.183 0.0100 
0.8 0.2096 0.0088 0.559 0.0049 1.205 0.0100 
0.9 0.2101 0.0093 0.553 0.0070 1.213 0.0099 
1.0 0.2090 0.0095 0.558 0.0072 1.227 0.0098 



able X(t): 

X(t)~t A (l + A x r"/*), (19) 

where uj is a well-known exponent of corrections to scaling. This expression reflects the 
scaling transformation in the critical range of time-dependent corrections to scaling in the 
form of t~ w l z to the usual form of corrections to scaling t^ v in equilibrium state for time 



t comparable with the order parameter relaxation time t r ~ £ z Q(k£) 17L Field-theoretic 
estimate of the uj value gives uj ~ 0.80 in the two-loop approximation 15|. Monte Carlo 
study of Ballesteros and Parisi [51] shows that uj ~ 1.0. 

To analyze our simulation date we have used the linear approximation of the (Xt~ A ) 
on t~ w l z with the changing values of the exponent A and the exponent uj from the interval 
[0.7,1.0]. Then, we have investigated the dependence of the mean square errors a of this 
fitting procedure for the function Xt~ A (t~ ul / z ) on the changing A and uj. In Fig. H]we plot 
the a for the magnetization (Fig. HH), logarifmic derivative of the magnetization (Fig. Hb), 
and the cumulant (Fig. Hb) as a function of the exponents j3 /uz, 1/vz, and d/z for uj = 0.8. 
Minimum of a determines the exponents z, u, and (3 for every uj. In Table U we present the 
computed values of the exponents (3/vz, l/vz, and d/z, and minimal values of the mean 
square errors a in these fits for values of the exponent uj = 0.7, 0.8, 0.9, 1.0. We see that 
the values of /3/vz, l/vz, and d/z are weakly dependent on the change of the exponent uj 
in the interval [0.7,1.0], but the uj = 0.8 is preferable because it gives the best fit for the 
magnetization and the logarifmic derivative of the magnetization dates. Finally, for the 



uj = 0.8 we find the following values of the exponents 



z = 2.489 ±0.021, 

v = 0.719 ±0.022, (20) 
P = 0.375 ± 0.045. 



It is interesting to compare these values of exponents with those obtained in 13j with 
the use of the field-theoretic approach 

z = 2.495, 

v = 0.716, (21) 
P = 0.350, 

which demonstrate a very good agreement with each other, but show an essential difference 
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from Monte Carlo results of Ballesteros and Parisi [5] with v = 1.009(13) and (3 = 0.526(15). 

B. Evolution from a disordered state 

We have also performed simulations of evolution of the system with linear defects on the 
largest lattice with L = 128, starting from a disordered state with small initial magnetiza- 
tions m = 0.02 and m = 0.001 at the critical point. The initial magnetization has been 
prepared by flipping in an ordered state a definite number of spins at randomly chosen sites 
in order to get the desired small value of mo- In accordance with Section [Til a generalized 
dynamic scaling predicts in this case a power law evolution for the magnetization M(t), the 
second moment M^ 2 \t) and the auto-correlation A(t) in the short-dynamic region. 

In Fig. El and [7] we show the obtained curves for M{t) (Fig. EJ), M^(t) (Fig. [6} 
and A(t) (Fig. [7]), which are plotted in log-log scale up to t = 700. These curves were 
resulted by averaging over 3000 different samples with 25 runs for each sample. From 
Fig. [5] we can see an initial increase of the magnetization, which is a very prominent 



phenomenon in the short-time critical dynamics [17j, Il8| . But in contrast to dynam 



ics of the pure systems |l7j], we can observe the crossover from dynamics of the pure 
system on early times of the magnetization evolution up to t = 70 to dynamics of 
the disordered system with the influence of linear defects in the time interval [100,650]. 
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FIG. 5: Time evolution of the magnetization M(t) for L = 128 with the initial magnetization 
mo = 0.02 (a) and mo = 0.001 (b) at the critical temperature T c = 3.9281. 
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FIG. 6: Time evolution of the second moment for L = 128 with the initial magnetization 

mo = 0.02 (a) and mo = 0.001 (b) at the critical temperature T c = 3.9281. 



The same crossover phenomena were observed in 
evolution of the second moment (t) and the au- 
tocorrelation A(t). In result of linear approxima- 
tion of these curves in the both time intervals we 
obtained the values of the exponents 8, c<i and c a 
in accordance with relations in (jHJ), (GJ) and (ED for 
initial states with itlq = 0.02 and mo = 0.001 (Ta- 
ble HI]) • The final values of these exponents and also 
the critical exponents z, f3/v and x were obtained 
by extrapolation to m = 0. In Table HT1 we compare 
the values of these exponents with values of corre- 
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FIG. 7: Time evolution of the autocor- 
relation A(t) for L = 128 with different 
initial magnetization values mo at the 
critical temperature T c = 3.9281. 



TABLE II: Values of the critical exponents obtained in present work for evolution from disordered 
initial states with different too and extrapolated to mo = and corresponding exponents for pure 
systems \v\ and from fl^ 



e 


C2 


C a 


z 




Xq 


t € [10, 70] 



m = 0,02 0.086(12) 0.964(28) 1.384(26) 
to = 0,001 0.099(9) 0.973(19) 1.364(23) 

to = 0.101(10) 0.975(23) 1.363(26) 2.049(27) 0.501(27) 0.708(34) 



t € [100,650] 



to = 0,02 0.152(12) 0.812(21) 1.103(16) 

too = 0,001 0.149(10) 0.804(19) 1.047(12) 

too = 0.149(11) 0.801(20) 1.043(14) 2.517(32) 0.492(28) 0.867(37) 

TFD [12] 2.495 0.489 

pure [17] 0.108(2) 0.970(11) 1.362(19) 2.041(18) 0.510(14) 0.730(25) 



sponding exponents for the pure Ising model [18] and theoretical field description (TFD) 



results for system with linear defects [13J. The obtained values quite well agree with re- 



sults of simulation from an ordered state with mo = 1 and with results from 13( and show 
that LR-correlated defects lead to faster increasing of the magnetization in the short-time 
dynamic regime in compare with the pure system. 



C. Measurements of the critical characteristics in equilibrium state 

With the aim to verify the short-time dynamics method and the results obtained we also 
carried out the study of the critical behavior of 3D Ising model with the linear defects of 
random orientation by traditional Monte Carlo simulation methods in equilibrium state. 
For simulations we have used the Wolf single-cluster algorithm. We have computed for the 
critical temperature T c = 3.9281(1) the values of different thermodynamic and correlation 
functions in equilibrium state such as the magnetization, susceptibility, correlation length, 
heat capacity, and Binder cumulant U4 for lattices with sizes L from 16 to 128 for the 
same spin concentration p = 0.8. The use of well-known scaling critical dependences for 
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FIG. 8: Dependence of the mean square errors a of the fits for heat capacity (a), magnetization 
(b), susceptibility (c), and thermal derivative of cumulant (d) as a function of the exponents a/v, 
(3/v, j/u, and v for different values of u. 



these thermodynamic and correlation functions with taking into consideration the finite size 
scaling corrections 



C(L) ~ L a/u (1 + aL-") 
M(L) ~ L' p/v (1 + WT^) 
X (L) ~ (l + cL"") 



(22) 
(23) 
(24) 
(25) 



makes it possible to determine the critical exponents a, u, (3, 7, and to by means of sta- 
tistical data processing of simulation results. To analyze simulation data we have used the 
linear approximation of the (XL~ A ) on L~ u and then investigated the dependence of the 
mean square errors a of this fitting procedure for the function XL~ A (L~ UJ ) on the changing 



TABLE III: Values of the exponents a/v, (3/v, 
7/V, and v with values of the exponent uj, giv- 
ing the best fit in approximation procedure 



TABLE IV: Values of the critical exponents obtained 
in present work for average value of exponent uj = 
0.76 and corresponding exponents from 13J] 



a/v 0/v 7/V v a (3 7 v 

-0.096(3) 0.457(2) 2.032(1) 0.710(10) present -0.078(30) 0.362(20) 1.441(15) 0.710(10) 

uj 0.90 0.65 0.70 0.80 [13] -0.1048 0.3504 1.4453 0.7155 

exponent A and u values. In Fig. [8] we plot the a for heat capacity (Fig. Ek), magnetiza- 
tion (Fig. [8b), susceptibility (Fig. Eb), and temperature derivative of cumulant (Fig. EH) as 
a function of the exponents a/u, /3/V, 7/V, and v for different values of uj. Minimum of 
a determines the values of exponents. In Table II I II we present the obtained values of the 
exponents a/u, f3/v, 7/V, v, and uj, which give minimal values of a in these fits. Then 
we determine the average value of uj = 0.76(5) with the use of which there were computed 
the final values of exponents. In Table IIVI there are presented the values of the exponents 
obtained in this work by simulation methods and from [l^] with the use of the field-theoretic 
approach and scaling relations for critical exponents. 

The comparison these values shows their good agreement within the limits of statistical 
errors of simulation and numerical approximations and good agreement with the values of 
the static critical exponents computed by the short-time dynamics method. 



IV. MEASUREMENTS OF THE CRITICAL TEMPERATURE AND CRITICAL 

EXPONENTS FOR 3D XY-MODEL 



Also, we have carried out the Monte Carlo 
study of the effect of LR-correlated quenched de- 
fects on the critical behavior of 3D XY-model 
characterized by the two-component order param- 
eter. As is well-known {2I, [3], renormalization 
group analysis predicts the possibility of a new 
type of the critical behavior for this model dif- 
ferent from the critical behavior of pure XY-like 
systems or systems with point-like uncorrelated 
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FIG. 9: Binder cumulant U 4 (T,L) for 3D 
XY-model ELS cl function of T for lattices 
with different sizes L. 



defects. We considered the same site-diluted cubic lattices with linear defects of random 
orientation in the samples with the spin concentration p = 0.8. The critical temperature 
T c = 1.8626(5) was determined by the calculation of Binder cumulant Ui(L,T) for lattices 
with sizes L from 32 to 128 (Fig. [9]). For simulations we have used the Wolff single-cluster 
algorithm. 



A. Evolution from an ordered state 



We have performed simulations of the criti- 
cal relaxation of the XY-model with linear de- 
fects starting from an ordered initial state. As 
example, in Fig. [TU] we show the obtained curves 
for the magnetization M(t) (Fig. [TUh). Binder 
cumulant Ui(t) (Fig. HDb) and the logarithmic 
derivative of the magnetization d T lnM(t,r)\ T=0 
(Fig. [TUb). which are plotted in log- log scale up 
to t = 1000 for lattices with L = 128. These 
curves were resulted by averaging over 3000 dif- 
ferent samples. On these figures we can observe 
the crossover from dynamics, which is similar to 
that in the pure system on early times of the evo- 
lution up to t = 150, to dynamics of the dis- 
ordered system with the influence of linear de- 
fects in the time interval [350,800]. The slope 
of M(t), U 2 (t) and d T lnM(t, r)| T=0 over the in- 
terval [350,800] provides the exponents /3/uz = 
0.239(1), d/z = 1.221(21) and 1/uz = 0.531(13), 
whereas our theoretical-field predictions [l3| give 
the following values of exponents fi/vz — 0.204, 
d/z = 1.268 and 1/uz = 0.556. As well as for Ising 
model, we have considered the corrections to the 
scaling in order to obtain accurate values of the 
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FIG. 10: Time evolution of the magne- 
tization M(t) (a), Binder cumulant U%(t) 
(b) and the logarithmic derivative of the 
magnetization d T lnM(t, t)\ t=q (c) for 3D 
XY-model with lattice size L = 128 at the 
critical temperature T c = 1.8626. 



critical exponents in concordance with procedure, which was discussed in Subsection IHI Al 
As a result of this analysis we obtained the following values of critical exponents 

z = 2.364 ±0.007, 
v = 0.778 ±0.026, 

(26) 

P = 0.370 ±0.030, 
to = 1.05 ±0.04. 

The comparison of these values of exponents with those obtained in [l^l with the use of 
the field-theoretic approach z = 2.365, u = 0.760, (3 = 0.366, and on = 1.15 in 2] shows 
their good agreement within the limits of statistical errors of simulation and numerical 
approximations. 

The obtained results confirm the strong effect of LR-correlated quenched defects on 
both the critical behavior of 3D Ising model and the systems characterized by the many- 
component order parameter. 



B. Evolution from a disordered state 



We have also performed simulations of evolution of the system with linear defects on 
the largest lattice with L = 128, starting from a disordered state with small initial mag- 
netizations m = 0.01, m = 0.0075 and m = 0.005 at the critical point. The initial 
magnetization has been prepared by flipping in an ordered state a definite number of spins 
at randomly chosen sites in order to get the desired small value of m . In accordance with 
Section [Til a generalized dynamic scaling predicts in this case a power law evolution for 
the magnetization M(i), the second moment M^ 2 \t) and the auto-correlation A{t) in the 
short-dynamic region. 

In Fig. [11] and [12] we show the obtained curves for M(t) (Fig. HT), M< 2 )(t) and A(t) 
(Fig. IT2"]) . which are plotted in log-log scale up to t = 700. These curves were resulted by 
averaging over 3000 different samples with 25 runs for each sample. From Fig. [TT] we can 
see also an initial increase of the magnetization, which is a very prominent phenomenon 
in the short-time critical dynamics. But in contrast to dynamics of the pure systems, we 
can observe, as previously for Ising model, the crossover from dynamics of the pure system 
on early times of the magnetization evolution up to t = 100 to dynamics of the disordered 
system with the influence of linear defects in the time interval [200,650]. In this time interval 
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FIG. 11: Time evolution of the magnetization M(t) for 3D XY-model with the initial magnetization 
mo = 0.01 (a), mo = 0.0075 (b) and mo = 0.005 (c) at the critical temperature T c = 1.8626; 
determination of asymptotic value of the dynamic critical exponent 8 in the limit mo — > (d). 

we determined the values of the dynamic critical exponent 9, which are equal 9 = 0.182(3) 
for case with the initial magnetization mo = 0.01, 9 = 0.232(7) for m = 0.0075 and 
9 = 0.278(11) for m = 0.005. Then, in accordance with the asymptotic value of the 
exponent 9 = 0.374(14) was determined in the limit m — > (Fig. [TTH) with the use of linear 
extrapolation. 

The analysis of results for evolution of the second moment M^\t) (Fig. H~2k) and the 
autocorrelation A(t) (Fig. [I2"b). obtained for simulation of systems with the initial magne- 
tization m = (in reality for m = 10 -6 such as for XY-model the spin configuration with 
m = is impossible to prepare), gives directly the values of exponents c 2 = 0.825(23) and 
c a = 0.907(30) for the time interval [300,650]. The same crossover phenomena is observed in 
evolution of M^ 2 \t) and A(t) from dynamics of the pure system on early times to dynamics 
of the disordered system with the influence of linear defects. 

On basis of these values of exponents 9, c 2 and c a we obtained the exponents z = 2.342(57) 




FIG. 12: Time evolution of the second moment (a) and the auto-correlation A(t) (b) for 3D 
XY-model with the initial magnetization niQ = at the critical temperature T c = 1.8626. 

and (3/v = 0.534(35), which quite well agree with results of simulation from an ordered 
initial state with m Q = 1 and with results of theoretical field description for XY-like systems 
with linear defects [131] within the limits of statistical errors of simulation and numerical 
approximat ions . 



V. CONCLUSION REMARKS 



The present results of Monte Carlo investigations allow us to recognize that the short- 
time dynamics method is reliable for the study of the critical behavior of the systems with 
quenched disorder and is the alternative to traditional Monte Carlo methods. But in contrast 
to studies of the critical behavior of the pure systems by the short-time dynamics method, 
in case of the systems with quenched disorder corresponding to randomly distributed linear 
defects after the microscopic time t m i C ~ 10 there exist three stages of dynamic evolution. 
For systems starting from the ordered initial states (mo = 1) in the time interval of 50-200 
MCS, the power-law dependences are observed in the critical point for the magnetization 
M(t), the logarithmic derivative of the magnetization d T lnM(t, t)| and Binder cumulant 
U2(t), which are similar to that in the pure system. In the time interval [450,900], the power- 
law dependences are observed in the critical point which are determined by the influence 
of disorder. However, careful analysis of the slopes for M(t), d T lnM(t,r)\ T=0 and C/ 2 (t) 
reveals that a correction to scaling should be considered in order to obtain accurate results. 
The dynamic and static critical exponents were computed with the use of the corrections 



to scaling for the Ising and XY models with linear defects, which demonstrate their good 
agreement with results of the field-theoretic description of the critical behavior of these 
models with long-range correlated disorder. In intermediate time interval of 200-400 MCS 
the dynamic crossover behavior is observed from the critical behavior typical for the pure 
systems to behavior determined by the influence of disorder. 

The investigation of the critical behavior of the Ising model with extended defects starting 
from the disordered initial states with mo — also have revealed three stages of the dynamic 
evolution. It was shown that the power-law dependences for the magnetization M(t), the 
second moment M^it) and the autocorrelation A{t) are observed in the critical point, which 
are typical for the pure system in the interval [10,70] and for the disordered system in the 
interval [100,650]. In intermediate time interval the crossover behavior is observed in the 
dynamic evolution of the system. The obtained values of exponents demonstrate a good 
agreement within the limits of statistical errors of simulation and numerical approximations 
with results of simulation of the pure Ising model by the short-time dynamics method [l^] 
for the first time interval and with our results of simulation of the critical relaxation of this 
model from the ordered initial states. 

Also, we would like to note that over complicated critical dynamics of the systems with 
quenched disorder the accurate determination of the critical temperature is better to carry 
out in equilibrium state from the coordinate of the points of intersections of the curves spec- 
ifying the temperature dependence of Binder cumulants U^L, T) or ratio £/L for different 
linear sizes L of lattices. 

The obtained results for 3D XY-model confirm the strong influence of LR-correlated 
quenched defects on the critical behavior of the systems described by the many-component 
order parameter. As a result, a wider class of disordered systems, not only the three- 
dimensional Ising model, can be characterized by a new type of critical behavior induced by 
quenched disorder. 

We are planning to continue the Monte Carlo study of critical behavior of the model with 
LR-disorder for different spin concentrations p and investigate the universality of critical be- 
havior of diluted systems with LR-disorder focusing on the problem of disorder independence 
of asymptotic characteristics. 
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